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Abstract 
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new mass parameter in the heavy quark propagator. 
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1 Introduction 



The Isgur-Wise (IW) function £(o>) [T], the universal function describing a class 
of B(vb) — > D(v c ) transitions has been studied for many years. (Here Vf> and v c 
are the four-velocities for the mesons containing a b, or a c-quark, respectively, 
and uj = Vb ■ v c ). While some general features of the function follow from heavy 
quark symmetry [2], its more detailed shape has been studied in various model 
approaches [SUlElElIZlElEKinKII]. 

It is well known that understanding the shape of the IW function, and in 
particular its slope £'(1) at the zero recoil point uj = 1 is a necessary prerequisite 
for determination of the V c b element of the CKM quark mixing matrix. However, 
there is also another incentive to study this function in quark models — namely, 
it is an essential ingredient in the model description of amplitudes for important 
non-leptonic decays of the type B — > DD, DK, Dir. Thus, to approach such 
theoretically difficult processes, one first needs a reasonable description of the IW 
function. 

Some versions of quark models [31 HI [7] gave a slope £'(1) — — 0.4 to —0.3, which 
is not in agreement with general theoretical expectations expressed in Bjorken 
[j~2] and Uraltsev [13] sum rules which together imply —£'(1) — 3/4. Also, a 
combined fit [T4] to results of experimental measurements of B — > D*lv decays 
gives —£'(1) = 1.16 ± 0.05, and, although dispersion of experimental results is 
large leading to a small confidence level of this fit (« 1%), it seems reasonable 
to assume that absolute value of the slope cannot be significantly smaller than 1. 
In this paper we propose a modified version of the model in [TJ, which has a 
particular feature of explicit inclusion of the gluon condensate effects, enabling 
consistent estimation of non-factorizable amplitudes [T51 [TBI [UJ (TSJ [T9], and we 
demonstrate that this model gives a satisfactory description of the IW function 
slope. 

There are two slightly different philosophies within the family of heavy-light 
chiral quark models. In Refs. [31 SI EJ the focus is on the bosonization procedure. 
The quark Lagrangian is bosonized by attaching meson fields to quark loops at 
zero external momentum, thereby integrating out the quarks. External momenta 
then correspond to derivatives of meson fields. On the other hand, in the approach 
of the Bari group [6] the external momenta are kept in the loop integral and mass 
differences between heavy mesons and the heavy quark appear in the final result. 
This mostly works fine in [6j, but a problem seems to arise when one tries to 
describe transitions between heavy mesons with different masses. 

In this paper we will work with zero external momenta (focus on bosonization), 
but introduce a parameter A in the heavy quark propagator, which will, for positive 
values of A, correspond to an extra dynamical mass of the heavy quark. 

In Sect. H]we describe this extended version of heavy-light chiral quark model. 
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In Sect. |3] we show how the model is bosonised and how the model parameters 
are related to physical quantities. In Sect. H] we calculate the IW function, and 
then we conclude. Appendix |A] contains recursion formulae and expressions for 
the relevant heavy-light loop integrals. 

2 Heavy-Light Chiral Quark Model (HL^QM) 

The total Lagrangian describing both quark and meson fields is [7j: 

£ = C-HQET + £-xQM + C-Int , (1) 

where [2] 

Chqet = ~Q v {iv-D - A) Q v + 0(m Q l ) (2) 

is the Lagrangian for Heavy Quark Effective Field Theory (HQEFT), with the 
mentioned extra mass added. The heavy quark field Q v annihilates a heavy quark 
with velocity v and mass vhq. Moreover, is the covariant derivative containing 
the gluon field (eventually also the photon field). In the C^m^ 1 ) term was 
also considered, but it will not be needed in this paper, because l/m c and l/m b 
corrections will not be calculated. 

The light quark sector is described by the Chiral Quark Model (xQM), having 
a standard QCD term and a term describing interactions between quarks and 
pseudo-scalar light mesons: 

£>xQM = q{il^D^ - M q )q - m(q R T,^q L + q L T,q R ) , (3) 

where q T = (u, d, s) is the light quark field triplet. The left- and right-handed 
projections and are transforming after SU(3)l and SU(3)r, respectively. 
M. q = diag(m u , m^, m s ) is the current quark mass matrix, m is the (577(3) invari- 
ant) dynamical mass of light quarks, and £ = exp(2«n// 7r ), where IT is a 3 by 3 
matrix containing the pseudo-scalar meson octet (jr, K, rj) in a standard way. 

There is also a "rotated version" of the xQM with flavour-rotated quark fields 
X given by: 

xl = ; xr = €qr ; f • f = s • ( 4 ) 

(This field £ containing the light mesons should be distinguished from the IW 
function In the rotated version, the chiral interactions are transformed into 

the kinetic term, while the interaction term proportional to m in Q becomes a 
pure (constituent) mass term [20| [T5]: 

£xQM = X [lUiDp + + 7s A) " m \ X ~ xM qX , (5) 
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where the vector and axial vector fields and A a are given by: 



and 

M q = ?M\£L + ZM\£R. (6) 

Here L is the left-handed projector in Dirac space, L = (1 — 7s)/2, and i? is the 
corresponding right-handed projector. The Lagrangian ([5]) is manifestly invariant 
under the unbroken SU(3)v symmetry. In the light sector, the various pieces 
of the Lagrangian describing strong interactions of mesons can be obtained by 
integrating out the constituent quark fields \i an d these pieces can be written in 
terms of the fields A^ and M q which are manifestly invariant under local SU(3)v 
t r ansf or mat ions . 

In the heavy-light case, the generalization of the meson-quark interactions of 
the pure light sector xQM is given by the following SU(3) invariant Lagrangian: 

Hint = ~G H [XaHSQv +Q~ V H a vX a] + \H* H$\ , (7) 

where Gh and G3 are coupling constants, and H" is the heavy meson field con- 
taining a spin zero and spin one boson: 

H a v = p + {p;r - *p 5 a 7 5 ) , 

HS = 7°(^ a ) t 7° = [{P;) ] l» - i(P 5 a H] P + , (8) 

where 

P± = 1 ± 2 7 ' V ■ H vl -v=-H v , 7 -vH~ v = -ir v . (9) 

The field P^{P^) annihilates a heavy-light meson, 0~(1 - ), with velocity v. The 
index a runs over the light quark flavours u, d, s, and the projection operators P± 
have the property 

P±l fM P± = ± P± P± . (10) 

Note that in Refs. [21 HI El [5] Gh = 1 is used. However, in that case one uses a 
renormalization factor for the heavy meson fields H v , which is equivalent to the 
approach in [7] and here. The term oc I/G3 is (partially) cancelled by a self-energy 
loop of order G H , and determines the mass difference between the heavy quark 
and the corresponding heavy meson. 

In our model, the hard gluons are considered to be integrated out and we are 
left with soft gluonic degrees of freedom. These soft gluons can be described us- 
ing the external field technique, and their effect will be parameterized by vacuum 
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Figure 1: Feynman rule for the light quark - soft gluon vertex. 



expectation values, such as the gluon condensate (^G 2 ). Gluon condensates with 
higher dimensions could also be included, but we truncate the expansion by keep- 
ing only the one with lowest dimension. When calculating the soft gluon effects in 
terms of the gluon condensate, we follow the prescription given in [21]. The cal- 
culation is easily carried out in the Fock-Schwinger gauge, where one can expand 
the gluon field as 

A a »(k) = -^G^(o) AjW^) + " " • • (11) 

Since each vertex in a Feynman diagram is accompanied by an integration, we get 
the Feynman rule given in figure [TJ The gluon condensate is obtained by averaging 
in colour space which yields the following replacement rule: 

9 s 2 GpG^ -> — - - ^J) ab (-^-G 2 )—{g lia g v p - g^g va ) ■ (12) 



3 Bosonization within the HLxQM 

The interaction term Cj nt in (J7]) can now be used to bosonize the model, i.e. to 
integrate out the quark fields. This can be done in the path integral formalism and 
the result is formally a functional determinant. This determinant can be expanded 
in terms of Feynman diagrams, by attaching the external fields H®, H£, V M , and 
M q of section E] to quark loops. Some of the loop integrals will be divergent and, 
analogously to the pure light sector case [201 E21 [221 EES] , they have to be related 
to physical parameters. The resulting strong chiral Lagrangian has the following 

form [22 EH ESI EH EH E3 ED]: 

Cstr = -Tr \H~ a {iv ■ V)H a ] - g A Tr [~H~ a H bl ^ 5 Al] + .... , (13) 

where the dots indicate other terms of higher order in the chiral expansion, and 
the covariant derivative T> contains both the photon field and the field V. The 
1 / rriQ suppressed terms have been discarded in the present paper. 
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Va 5 Aa 



Vq , Aa 




Figure 2: Coupling to vector and axial vector current 



The Feynman diagrams responsible for the kinetic and axial vector terms in 
(1131) are shown in Fig. [2J As mentioned in the Introduction, these two diagrams 
are calculated at zero external heavy meson momentum. The non-gluonic loop 
integral (first on Fig. [2]) for the strong vector or axial vector current is 

J x = -N C J ^Tr{(-iG H W v )iS v (k)(-iG H H v )iS(k)V x iS(k)] , (14) 

where Ty = 7 M and = 7^75 for couplings to X = V and X = A. Further, 
S v (k) and S(k) are the heavy quark propagator and the standard light quark 
propagator, respectively: 

Sv{k) = (Vk-A) = (fc2 + + m 2) • ( 15 ) 

In previous papers [S HI [7J [6] , A = was assumed, but here we let A^O. 

The gluonic part of the bosonized currents for one of the diagrams (lower left 



6 



on Fig. ED is 

J XjGl = - J -^Tr[(-iG H lT v )iS v (k)(-iG H H v )iS(k)V x iS(k) 

[ig s l a A a {q 2 )\ iS{k - q 2 ) [ig sl a A a { qi )} iS{k - q x - g 2 )} , (16) 



where gluon fields are represented by the expression from Eq. flTTj) . and it is un- 
derstood that the derivatives with respect to soft gluon momenta are to be applied 
to the whole integrand. There are two more diagrams with different ordering of 
gluon and (axial) vector vertices and after adding all four diagrams from Fig. [2] we 
obtain 



^,Tot = -9xTr{H v H v T^} , (17) 



where 



~ 2 ■■ "> v 



9x = iG H N c <i R x - — Z X (-^G') 



R v = -2(m - A)I 2 - I 1A - 2A(m - A)J 2 ,i , (19) 
Z v = UAmh + 192m 2 (m - A)J 5 + 24m(m + 6A)/ 4 ,i 

+ 192m 2 A(m - A)/ 5i i , (20) 
and for the axial case 

Ra = -^(3m - A)/ 2 + -I ltl + H( m - A) (2m - A)J 2>1 , (21) 

Z A = A8ml 4 + 64m 2 (3m - A)J 5 - 8m(13m - 6A)/ 4 ,i 

-64m 2 (m- A) (2m- A)J 5)1 . (22) 

The loop integrals /„ and 7 nr occurring in the expressions above are defined in 
Appendix |Al where also expressions for the finite ones are given. The integrals 7 2 , 
Ix t i above and I\ from Eq. (|43|) below are logarithmically, linearly, and quadrati- 
cally divergent, respectively. They will be expressed in terms of model parameters. 
The negative parity axial coupling constant g A is taken as model input parameter, 
9a = 0.59, and gy = 1 (i.e. the normalization). 

Within the pure light sector the logarithmically and quadratically divergent 
integrals are related to the pion decay constant f n and the quark condensate ( qq ) 
in the following way [221 [231 [15] : 

f^-iim'NJz + ^-^C?), (23) 
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1 n 

{qq) = -AimNJ, - — — (— G 2 ) . (24) 

12m 7i 

This is obtained by relating loop diagrams to physical quantities analogously to 
Eqs. (1291) and (|30|) below. (Here the a priori divergent integrals I± and Ii have to 
be interpreted as the regularized ones.) Since the pure light sector is a part of our 
model, we keep these relations in the heavy-light case studied here. In addition, 
in the heavy-light sector the (formally) linearly divergent integral I\ i will also 
appear. It will be related to the physical value of using Eq. ([TBI) for X = A. 

Eliminating thus I\ y \ from the ([TBI) and inserting the expression for I2 obtained 
from f[23[) we find the following expression for Gh'- 

2m . . 

G 2 H = — PA , (25) 

J 7T 

where the quantity pa is of order one and given by 

Pa = — 7 x r > I 26 ) 

4 (1 - & + - a^c 2 ) 



where 



and 



k a = il67rm(l - — ) 2 h 1 , (27) 
m 



A 1 

* = 12 



(2J 



Let us now consider these relations in two characteristic limiting cases: A — > 
and A — > m. 

In the limit A — > 0, f[T8j) reduces to 



where 



1 = -iG 2 H iV c {/ 3/2 + 2mJ 2 + ^(-G 2 )} , (29) 
,^^{1^-2^-^-^^)}, ,30) 

^3/2 = (A,i)a->o ; «v = - «U = • (31) 

This is the result of [7j, except for the sign of ka which is wrong there. Numerically, 
this change has no dramatic consequences. In the limit A — > we also have 

ka -> 1 ; V? ^ Vi = j^- i^v + 3k a ) = • (32) 
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Figure 3: The decay constant /# in dependence on m and A. The condensates 
are taken to be (qq) = (-0.27 GeV) 3 and (^G 2 ) = (0.32 GeV) 4 . 



Eliminating I2 we obtain the relation for I3/2 
- iNJ 3/2 



3(1 -fiu) m 3(k a - k v ) ,^s_q2\ 



16tt 



4m 3 



7T 



(33) 



which will replace the Eq. (39) of [?]. 

In the limit A — > m, which means that heavy and light quarks have the same 
constituent mass, we obtain 



9A = G%N Jl iNJu -l ml2 -^ { ^ )} . 



(34) 
(35) 



Then, for A — > m, there is a simplification because ka — > and 77^ — ^ and 
thereby 

PA^\(l + Sg A ). (36) 

The gluon condensate may be related to the matrix element of the chromo- 
magnetic interaction [2]: 

3A 2 = &{H) = ^M(H\Q v ^a • GQ V \H) = \m Q {M H * - M H ) . (37) 

Such a link was used in [7], but we will not use it because it formally belongs to 
1/rriQ corrections, which are not considered here. Also, it turns out that such a 
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Figure 4: Diagrams for bosonization of the left handed quark current to 
leading order, determining a# 

choice doesn't lead to any dramatic differences, when numerical results in A — > 
limit are compared to those from [7J. 

Within the full theory (Standard Model) at quark level, the weak current is 

J? = q fl a (l- l5 )Q, (38) 

where Q is the heavy quark field in the full theory. Within HQEFT this current 
will, below the renormalization scale fi = mg (= m&,m c ), be modified in the 
following way [2]: 

J? = Xa£ f F a Qv + 0(m Q l ) , (39) 

where 

T a = C~ tl a L + C v v a R. (40) 

Bosonising this current by calculating diagrams from Fig. H]we obtain to zeroth 
order in the axial field and to first order in the gluon condensate the weak current 



ja _ ^H_rp 

J f ~ 2 1 r 



(41) 



where heavy meson decay constant is given by 



h = vmf a " • (42) 

with C 7 = 1.077 and C v = 0.0489 being Wilson coefficients within HQEFT [2]. In 
addition /# and fo have chiral corrections of order 20 MeV (see [7] and references 
therein), as well as I/tuq corrections not considered here. From the diagrams in 
Fig. HI we obtain: 

a H = -2iG H N c { - I x + (m - A)/ a)1 

+ ^(-G 2 ) [-mh + J 3|1 + m(m - A)/ 4>1 ] } , (43) 
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where divergent integrals I\ and I^i will be expressed in terms of model parame- 
ters. Since this is the only predicted quantity depending on Ii it is easy to accom- 
modate physical values of fn, using (124)) and adjusting the quark condensate, see 
Fig. 1(b). In the limit A -»• 

a H = -2iG H N c (~h + mI 3/2 + l^L-^G 2 )] . (44) 
\ JS C m l 7r J 

where = (37r — 4)/384, in agreement with [7]. Furthermore, eliminating diver- 
gent integrals, one obtains 

Gh ( (qq) m 2 3(1 - g A ) 2 (3tt - 8)) a s 2 A 

In the limit A — > m 

a H ee -2iG H N c (-1, + ^^(-G 2 )) , (46) 
\ 96iV c m z tt J 

or, eliminating divergent integrals, 

This completes the specification and bosonisation of the HLxQM. The remain- 
ing free parameters of the model are two condensates (qq) and (^G 2 ) and two 
constituent masses m and A. We can now apply the model to calculation of 
phenomenological quantities, starting with the Isgur-Wise function. 



4 The Isgur-Wise function 

The Isgur-Wise function pQ, £(u), relates all the form factors describing the pro- 
cesses B — > D(D*) in the heavy quark limit. In our framework it can be defined 
via bosonisation of heavy-heavy quark current responsible for B — > D transition: 

= JLc ■ (48) 

Here Q Vc and Q Vc are the c and b quark fields within HQEFT. The IW function 
can be determined by calculating the diagrams shown in figure 

One normally expects that emission of the soft gluons from a heavy quark 
doesn't occur at the zeroth order in 1/mg. Namely, using the HQET Lagrangian 
(j2J), and differentiating the expressions involving heavy quark propagators accord- 
ing (|11|) . will naturally, for diagrams of same class as those on Fig. HJ lead to 



(+) 



-£{u)Tr 



(+) 
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(a) 






p — ki — k'2 




+ ki + k 2 



Figure 5: Bosonisation corresponding to the Isgur-Wise function. 



expressions proportional to v^G® = 0, where t> M is either Vb or v c . However, for 
the Isgur-Wise function there are two velocities (t>& and v c ) in play in the diagram. 
Therefore one may obtain contributions proportional to 



(49) 



which will generally, away from the strict heavy quark limit w — > 1, be different 
from zero. Let us also mention that some care is required because momentum 
flow in the diagrams for the translationally non-invariant amplitudes in the Fock- 
Schwinger gauge is non-trivial, see Fig. [SJ 

The corresponding results for diagrams (a)-(d) are 



£(w) a = iG 2 H N c (m - 



uj + 1 



A)/ 



1,1,1 



u + 1 



h. 



aco) b = zG 2 H (-G 2 )mrf 2 (m 2 h^ + J 3)ljl - (7 4 .i + A/^,^ 

7T \ Co" + 1 / 



^(^GV7r 2 (a;-l)J 2j2j2 
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G 2 H (^G 2 >7r 2 (a; - 1) ( - J li2)3 - Ii, 3 ,2 + (m(a> + 1) - 2A)/ W 



(50) 
(51) 
(52) 

(53) 
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Figure 6: The slope of Isgur-Wise function £'(1) in dependence on the gluon 
condensate for various choices of the light quark constituent mass m (a), 
and the decay constant fs in dependence on the quark condensate (b), for 
A = 0.1 GeV. Complementary plots are not so interesting because the 
dependence of £'(1) on the quark condensate and the dependence of fs on 
the gluon condensate are rather small. 

Since loop integrals have at the least two heavy and one light quark propagator, 
integrals of the type I n ,r,s occur above, and they are defined in the Appendix |A] 
We find that the identity = 1) = 1 follows from normalization of the vector 
current in eq. flTSJ): 

£(u = 1) = (£, + & + £ c + £ d ) w=1 = (£, + &) w=1 = 1 . (54) 

Finally, for the slope of the IW function in the no- recoil limit to — > 1, we have 

r(l)a = ^^ c Q/ lil -(im-^A)/ 1 , 2 -^A(m-A)(/ li 3-2m 2 / 3 , 1 )^ , (55) 

ai) b = i G 2 H (^G 2 )mn 2 Q/ 2>4 + j(m - A)/ 3>4 + y foi + AI 4 , 2 )) , (56) 

(57) 



,2 ,^s n2 \ m ^ _ L 



71 iZ 



£'(1), = -iG 2 H (—G 2 )^-[ h i5 - (m - A)/ li6 

7T O 



(5? 



All the I n>r integrals above can be evaluated using formulas from Appendix |Aj 
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Figure 7: The slope of Isgur-Wise function £'(1) at no recoil point as a 
function of A for three values of m (a), and as a function of m for three 
values of A (b). The condensates are taken to be (qq) = (—0.27 GeV) 3 and 
(^G 2 ) = (0.32 GeV) 4 . 

5 Results and discussion 

Numerical results are presented in Figs. Eland CD Fig. E] displays the slope £'(1) 
of IW function as a function of the gluon condensate and the heavy meson decay 
constant /# as a function of the quark condensate, while Fig. [7] shows slope £'(1) 
for some generic condensate values in dependence on the dynamical masses A and 
m. One observes that for reasonable intervals of massed slope is of the order of 
-1 and beyond, which is in agreement with values mentioned in the Introduction, 
as well as with values obtained within other theoretical frameworks [21] and on 
the lattice [32]. It should be noticed that values for £'(1) °f order or bigger than 
one prefer smaller constituent mass than in [7j. For easier overview of analytical 
results, as well as for their numerical checks, it is convenient to again investigate 
limit A — > m. After using constant values for integrals in this limit, as given in 
Appendix [A] one obtains simple expression 

c , n s (39A ~ 1) (l + gftp f m 2 N c 11 ,(Xara\\ (kq\ 

e (1) " — i XfT~ \-^r + siE^ { ~ G } j ■ m 

Numerical values in this limit are not unreasonable. 

In conclusion, we have presented an improved Heavy-Light Chiral Quark Model 
where introducing additional mass parameter in the heavy-quark propagator re- 



3 Constituent mass of quarks in the presented model is smaller than in other similar quark 
models due to the explicit inclusion of gluon condensate in the dynamics of the model. 
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suited in a flexible model capable of consistent description of heavy-meson decays, 
where we placed particular emphasis on a characterization of Isgur-Wise function. 
Further applications of the model, such as calculation of non-factorizable ampli- 
tudes for non-leptonic heavy-meson decays could now be attempted. As the slope 
of the IW function is steeper than the one used in, say, Ref. [18], the partial 
amplitudes for B — > DD depending on the IW function, might be overestimated 
there. This will then have consequences for the size of the overall amplitude. 
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A Loop integrals 

Three-, two-, and one-point loop integrals with one light-quark propagator occur- 
ring in the calculations are: 

rap- = f pV •■• 

n ' r ' s J (2tt) 4 G» 2 - m 2 + ie) n {p ■ v - A + ie) r (p-v'-A + ie) s ' 1 ; 

rap- = f PV ••• (A2) 

J (2vr) 4 (p 2 - m? + ie) n (p • v - A + ie) r ' 1 ' 

rap- = f P°P P --- , A o^ 

J (2vr) 4 {p 2 - m 2 + ie) n ' 1 ' 
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To evaluate such integrals one first reduces tensor to scalar ones using relations 
(oj = v • v')\ 



1 n,l 



1 

" 3 
1 

2 
1 

+ 3 



4_ u + m 2 / na - 4A/„ - 4A 2 / na 



(A4) 



(A5) 



/ n _ 1 _ ( m 2 _ 4A 2 )4 _ 2AJ n _ M - 2A(m 2 - 2A 2 )/ n>1 
3 3 

-7 n _i + (-m 2 - A 2 )/„ + A/ n _i,i 



(A6) 



r 

n,r,s 



+ A(m 2 - A 2 )/ nil {v a g^ + + v 7 ^ a/3 ) 

l r 

; ^ (^7 n)rjS _i — J„ ir _i iS )f a + (W n)r _i jS — I n ,r,s-i)v a , 

+ A(^-l)/ nir > a + t; to )l . (A7) 



(Reduction of one-point 1% tensor integrals is simple and well-known.) Now 
all scalar two-point integrals I n , s can be reduced to linear combinations of I n> i 
integrals using general recursion formula 



-An 



s 2 - 3s + 2 



[n + s- 3)/ n+M _2 + m 2 (n + l)I n+2>s - 2 



s 2 -3s + 2 



+ A " [> , (A8) 



valid for s > 2, whereas the special case s = 2 is 

4,2 = -2n(/ n+ i + A7 n+lil ) . 



(A9) 



For the calculation of the slope £'(1) of Isgur-Wise function, one additionally 
needs derivatives of three-point integrals at point u = 1. Integrals themselves at 
this point are trivially given by 



7 n ,r,s(l) 



l n,r+s 



The derivatives are given by 



dl n ^ s (u) 



dui 



rs 



2(n- 1) 



I n -i,r+s+2 for n > 1 , 



(A10) 



(All) 
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while for n = 1 we have special case: 



dl lr>s (uj) 



dtu 



rs 



r + s + 



t( 2 



r + s — 1 + e 

r + s 



2m 2 
r + s 



l 2,r+s) , 



) 



(A12) 

where dimensional regularization parameter e is relevant only for divergent case 
r = s = 1 where using lim e _^ e/i j2 = —i/(8n) one gets 



du 



--(/ li2 + 2A/ li 3-4Am 2 / 3il ) . 



(A13) 



Then again recursion formulas above above can be used to reduce everything to 
I n i integrals. These integrals can be explicitly evaluated and they read 



In,l 

a 2 (x) 
a 3 (x) 
a±{x) 
a 5 (x) 
a 6 (x) 



lQn 2 m 2n 3 
1 



a n {—j 
m 



where 



1 — x 2 
1 1 



z->l 



4(1 -x 2 ) 2 
1 1 
~24(1 -x 2 ) 3 
1 1 



(F(x) - 2x + 2x 3 ) 



x-+l 1 

" 3 



(3F(x) - lOx + 14x 3 - 4x 5 ) 



2 
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(l5.F(x) - 66x + 118x 3 - 68x 5 + 16x 7 ) 
1920713^ ( 105 ^) - 558x + 1210x3 



192(1 -x 2 ) 4 
1 1 



(A14) 

(A15) 
(A16) 
(A17) 
^i>l (A18) 



— ~ • (A19) 



-1052x 5 + 496x 7 - 96x 9 ) 
Here the function F{x) is [33] 

^(x) = Vx 2 - 1 + ie ( log(x - Vx 2 - 1 + ie) - log(x + Vx 2 - 1 + k) ) (A20) 



which gives the various a n (0). Furthermore, JF(x) 



Note that ^"(x) n, 

—2xF(l/x), where F(l/x) is function from Eq. (A2) of 

This reduction of integrals is easy to implement on a computer and correspond- 
ing Mathematica code is available from the authors upon request. 
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